ROTATIONAL CORRECTION ON THE MORSE POTENTIAL THROUGH THE 



PEKERIS APPROXIMATION AND NIKIFOROV-UVAROV METHOD 

O . Ciineyt Berkdemir * 

Department of Physics, Faculty of Science and Literature, Erciyes University, 38039, Kayseri, Turkey 
Oh! and 

\ The Abdus Salam International Center for Theoretical Physics, Trieste, Italy 

o 

Jiaguang Han 1 

' Shanghai Institute of Applied Physics, Chinese Academy of Science, 201800, Shanghai, China 



(N 
00 



and 



The Abdus Salam International Center for Theoretical Physics, Trieste, Italy 

o ' 
in 
o 

S3 



Abstract 

The Nikiforov-Uvarov method is employed to calculate the the Schrddinger equa- 
tion with a rotation Morse potential. The bound state energy eigenvalues and the 
corresponding eigenfunction are obtained. All of these calculation present an effective 
and clear method under a Pekeris approximation to solve a rotation Morse model. 
Meanwhile the results got here are in a good agreement with ones before. 
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Introduction 



As an empirical potential, the Morse potential has been one of most useful and convenient 
model. The Morse potential gives a excellent qualitative description of the interaction be- 
tween the two atoms in a diatomic molecule and also it is a reasonable qualitative description 
of the interaction close to the surface jT|. For this potential, the Schrodinger equation shall 
be solved for the angular momentum quantum number / is equal to zero. As we know the 
rotation energy of a molecule is much smaller than that of vibration, and therefore in a pure 
Morse potential model the rotation energy of a molecule has been omitted. However, in some 
case, it is needed to be included if one wants to obtain analytical or semianalytical solutions 
to the Schrodinger equation. Then we will get a rotation Morse potential. Our interest is 
how to get the solution of Schrodinger equation for the rotation Morse potential. 

In the previous reports, several approximations have been developed to find better ana- 
lytical formulas for the rotating Morse potential j2]. Some of them require the calculation 
of a state-dependent internuclear distance through the numerical solutions of transcendental 
equations. Recently, an alternative method known as the Nikiforov-Uvarov (NU) method is 
introduced for solving the Schrodinger equation. The application of the method is to solve 
Schrodinger equation with some well-known potentials and to solve Dirac, Klein-Gordon 
and Duffin-Kemmer-Petiau equation for a Coulomb type potential. In the present work, the 
energy and the corresponding eigenf unctions are calculated employing the NU method under 
a Pekeris approximation [Hj. 

2 Basic Equations of Nikiforov-Uvarov Method 

The NU method provides us an exact solution of non— relativistic Schrodinger equation for 
certain kind of potentials jlj. The method is based on the solutions of general second order 
linear differential equation with special orthogonal functions [5]. For a given real or complex 
potential, the Schrodinger equation in one dimension is reduced to a generalized equation of 
hypergeometric type with an appropriate s = s(x) coordinate transformation. Thus it can 
be written in the following form, 



where a(s) and a (s) are polynomials, at most second— degree, and r (s) is a first— degree 
polynomial. To find a particular solution of Eq.([T|) by separation of variables, we use the 
following the transformation 




(1) 




(2) 
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This reduces Schrodinger equation, Eq.(l), to an equation of hypergeometric type, 



a(s)y" + r(s)y' + \y = 0, (3) 

where 0(s) satisfies = ir(s)/a(s). y(s) is the hypergeometric type function whose 

polynomial solutions are given by Rodrigues relation 

^> = ^™ s >]> < 4 > 

where B n is a normalizing constant and the weight function p must satisfy the condition jlj 

(ap)' = Tp . (5) 

The function it and the parameter A required for this method are defined as 



7T = ± 



\ 



2 

- a +ka (6) 



and 

A = k + tt'. (7) 

Here, tt(s) is a polynomial with the parameter s and the determination of k is the essential 
point in the calculation of tt(s). Thus, in order to find the value of k, the expression under 
the square root must be square of a polynomial. Hence, a new eigenvalue equation for the 
Schrodinger equation becomes 

/ n(n — 1) „ 

X = \ n = -nr'- 1 2 V , (n = 0,l,2,...) (8) 

where 

r(s) =r (s) + 2tt(s), (9) 
and it will have a negative derivative. 

3 Rotational Correction on the Morse Potential 

Consider a diatomic molecule system with reduced mass p in the Morse potential 

V{r) = D[e- Mr ~ ro) - 2e- a{r ~ ro) ] (D > 0, a > 0), (10) 

where D is dissociation energy, r is the equilibrium bond length and a is a parameter 
controlling the width of the potential well. The vibrations of a two-atomic molecule can be 
excellently described by this potential type and solved bound states for / = 0. If one want 
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to obtain the solution of I 7^ 0, the centrifugal term has to be attached to the potential in 
the Schrodinger radial equation. In this case, total potential shape becomes as follows 

WO = D[e- 2a ^ - 2e~ a ^} + ^ \ l) , (11) 

2^r z 

defining as the sum of the Morse potential and the centrifugal barrier. In order to calculate 
the energy eigenvalues and the corresponding eigenfunction, the potential function given by 
Eq. lfTUj) is inserted into the radial Schrodinger equation 

( %% * + D[e-**-**) - 2e-«*-"»>] + ^ Rm(r) = E nl R nl (r). (12) 



\ 2\i dr 2 L 2fir 2 

An analytical exact solution of this differential equation cannot be found without an approx- 
imation. For this case, we now outline the procedure of Pekeris [3J E] . 

3.1 Overview of the Pekeris Approximation 

The approximation is based on the expansion of the centrifugal barrier in a series of ex- 
ponentials depending on the internuclear distance, until the second order. However, by 
construction, this approximation is valid only for lower vibrational energy states. Therefore, 
for a Pekeris approximation, we can take care of the rotational term in the following way. 
We let x = (r — r )/r and around x — it may be expanded into a series of powers as 

V r<*( x ) = 77^^ = 7(1 - 2x + 3x 2 - Ax 3 + ...), (13) 

(i + xy 



7 = ^^, (M) 



with 

h 2 1(1 + 1) 

the first few terms should be quite sufficient. Instead, we now replace the rotational term 
by the potential 

V rot (x) = 7 (Dq + D ie ~ ax + D 2 e~ 2ax ) , (15) 
where a = ar . Combining equal powers of Eqs.([13p and (|T5j) we get 



~ , . ,., a 2 x 2 a 3 x 3 . ,„ n Aa 2 x 2 8a 3 x 3 . 
V rot (x) = 1 lD + D 1 (l-ux + — — + ...) + D 2 (l - 2ax + — — + ...) 



(16) 



Vrot(x) = 1 \D Q + D 1 + D 2 - x{D ia + 2D 2 a) + ^Piy + 2,02a 2 ) 



x 3 (D 1 — + D 2 — 
6 3 



+ 



(17) 
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where 



D = l-- + - 2 

a cr 

4 6 

a a 2 
1 3 

— + — • 

a cr 



D 1 



;is) 



We now can take the potential V rot instead of the true rotational potential V rot and solve the 
Schrodinger equation for I ^ in Eq. (fI2j) . 

3.2 Applying Nikiforov-Uvarov Method to the Schrodinger Equation 

Now, in order to apply the NU— method, we rewrite Eq. ()12j1 by using a new variable of the 
form s = e~ ax , 



d 2 R nl (s) ldR nl (s) 2fir 



ds 2 



+ 



s ds 



+ 



h 2 a 2 s 2 



(E nl - 7 A)) + (2D - 7 £>i)s ~(D + 7 D 2 )s : 



Rrd(s) = 0. 



By introducing the following dimensional parameters 



2nrl(E nl - 7 £>c 



01 " h 2 a 2 ' £2 h 2 a 2 

which leads to a hypergeometric type equation defined in Eq.((TJ): 

d 2 R n i , 1 c^-Rni 1 r 2 , 

~~nr H 3 — h ~j x _e i + - / 

ds z s ds s z L 
After the comparison of Eq.(|21|) with Eq. ((!]), we obtain the corresponding polynomials as 



^3 



2^r 2 {D + 1 D 2 ) 
h 2 a 2 



(19) 



(20) 



Rnl — 0. 



(21) 



r fs) = 1, a(s) 



a (s) = -ej + e 2 s - e 3 s 2 . 



[22) 



Substituting these polynomials into Eq.©, we obtain n function as 



7T 



s) = ±Je 3 s 2 + (k - e 2 )s + e\ 



(23) 



taking <j'(s) = 1. The discriminant of the upper expression under the square root has to be 
zero. Hence, the expression becomes the square of a polynomial of first degree; 



(k - e 2 ) 2 - 4e 2 e 3 = 0. 



(24) 



When the required arrangements are done with respect to the constant k, its double roots 
are derived as k + - = e 2 ± 2e\^fe~l. 

Substituting k + _ into Eq.([23|). the following possible solutions are obtained for ir(s) 
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7T(s) = ± 



(25) 



(V^3 s + £ l)> for K = E 2 + 2e ly /s^ 

It is clearly seen that the energy eigenvalues are found with a comparison of Eq.(J7J) and 
Eq.fjHJ). From the four possible forms of the polynomial 7r(s) we select the one for which the 
function r(s) in Eq.Q has a negative derivative. Therefore, the function r(s) satisfies these 
requirements, with 



t(s) = 1 + 2ei - 2^ s, 
r'(s) = -2^ . 



From Eq.flBJ we also get 



and also 



A = e 2 - 2e\yfe- 



It is seen that the parameter —el has the following form 



- e. 



S2 (I 

— — — n H — 
L2v^ V 2 



(26) 

(27) 
(28) 

(29) 



Substituting the values of — ef, e 2 and e 3 into Eq.(j2HJ, we can immediately determine the 
energy eigenvalues E n i as 

1 2 

(30) 



h\l + 1) / 3 3 \ /i 2 a 2 

-^ni = ~ o I 1 1- 



where 



2^1 



So 



ar oPtq J 2/x 



— — — \ n H — 

2^ V 2 



2^ a 2 x/£s 



2/xD /(/ + 1) / 2 3 



2 1 9 2 



(31) 



The last equation indicates that we deal with a family of the rotating Morse potential. 
Of course, it is clear that by imposing appropriate changes in the parameters D, a and 
t~o, the rotational energy spectrum for the any molecules can be calculated by the Pekeris 
approximation and Nikiforov-Uvarov method as well as other methods |Hl El El • In this 
study, we calculate the rotating energy values of the Morse potential for the CO and LiH 
molecules. The explicit values of the energy for the different values of n and I are shown in 
Tables 1 and 2, for known values of their relevant potential parameters [TTl 1121 ITB*] . 

Let us now find the corresponding eigenfunctions for this potential. Due to the NU— method, 
the polynomial solutions of the hypergeometric function y(s) depend on the determination 
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of the weight function p(s) which is satisfies the differential equation [a(s)p(s)]' = r(s)p(s). 
Thus, p(s) is calculated as 

p (s) = s 1+2 ^e~ 2 ^ s . (32) 

Substituting into the Rodrigues relation given in Eq.(j3J), the eigenfunctions are obtained in 
the following form 

d n r 



y nl (s) = B nl s^ l+2 ^e 2 ^ 



,(n+l+2E!) -2y^J s 



(33) 



where B n \ is the normalization constant. The polynomial solutions of y n i{s) in Eq. (j3*3*j) are 
expressed in terms of the associated Laguerre Polynomials, which is one of the orthogonal 
polynomials, that is 

y nl (s) =4 +2ei ("), (34) 
where v = l^fel s. By substituting 7r(s) and a(s) into the expression </>(s)'/0(s) = it(s)/a(s) 
and solving the result differential equation, the other part of the wave function in Eq.(j2J) is 
found as 

<f>(s) = s ei e~^ s , (35) 

or in terms of v 

<f>{v) = (2 v ^)- £ V £l e-^ 2 . (36) 

Combining the Laguerre polynomials and 4>{v) in Eq.(J2J), the radial wave functions are 
constructed as 

R nl (r) = A nl (2^T 3 y^^e-»/ 2 Ll+ 2 ^(v), (37) 

where A n \ is a new normalization constant. It is clearly verified that R n i{r) satisfies the 
following requirement 

POO 

/ Rll( r ) dr < 0°- 
JO 

The constant A n \ is determined by making this integral equal 1, 



nl , 

a Jo 



Ll +2£ Hv) 



dv=\. (3* 



The integral in Eq. ()38|) can be evaluated by using the recursion relation for Laguerre poly- 
nomials and then the normalization constant can be found as 

, 2 _ 4cm!(l + n + £l ) 2 (2 v /£j) 2£ i 

Anl ~ (l + n + Jfej)! • {6J) 

Therefore, the simplest radial wave function becomes for n = 0: 



Mr) = 2 2 - £ Hl + e0y (1 " n 2£i) , e- u/2 - (40) 



4 Conclusions 



We have presented an approximation for the rotational correction on the Morse potential 
which leads to analytic calculations for the energy spectrums and wave functions. In these 
calculations, we have used a new method which is developed by Nikiforov-Uvarov and ap- 
plied the Pekeris approximation. Our main results are summarized in Eq. (fHU|) and Eq. (J37|) . 
This new method is tested by calculating the energies of some actual rotational states of 
the CO and LiH molecules and comparing the results with those of variational and shifted 
1/N expansion methods. The method used here is the best advantage than the other meth- 
ods, such as series solutions, shifted 1/N expansion, supersymmetric approach and Laplace 
transforms, due to the fact that a systematical one. 
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Table 1: Energy eigenvalues (in eV) for the different values of n and I for CO molecule, with 
D = 90540 cm' 1 , a = 2.2994 A -1 , r = 1.1283 A and p = 6.8606719 amu [TTllTSj . 



n £ NU Method Variational 1/N Expansion 



-11.091 -11.093 -11.091 

5 -11.084 -11.085 -11.084 

10 -11.065 -11.066 -11.065 



5 





-9.795 




-9.788 


5 


5 


-9.788 




-9.782 


5 


10 


-9.769 




-9.765 


7 





-9.299 




-9.286 


7 


5 


-9.292 




-9.281 


7 


10 


-9.274 




-9.265 
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Table 2: Energy eigenvalues (in eV) for the different values of n and I for LiH molecule, with 
D = 20287 cm- 1 , a = 1.1280 A -1 , r = 1.5956 A and n = 0.8801221 amu [TTllTSj . 



n t NU Method Variational 1/N Expansion 



-2.4287 -2.4291 -2.4278 

5 -2.4012 -2.4014 -2.3999 

10 -2.3287 -2.3287 -2.3261 



5 -1.6476 -1.6242 

5 5 -1.6236 -1.6074 

5 10 -1.5606 -1.5479 



7 -1.3774 -1.3424 

7 5 -1.3549 -1.3309 

7 10 -1.2957 -1.2781 
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